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Abstract 

Inclusive jet production, e + e~ — > e + e~ jet X, is studied using 560 pb _1 of data 
collected at LEP with the L3 detector at centre-of-mass energies between 189 and 
209 GeV. The inclusive differential cross section is measured using a k t jet algorithm 
as a function of the jet transverse momentum, p t , in the range 3 < p t < 50 GeV for 
a pseudorapidity, r], in the range — 1 < rj < 1. This cross section is well represented 
by a power law. For high p t , the measured cross section is significantly higher than 
the NLO QCD predictions, as already observed for inclusive and 7r° production. 
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1 Introduction 



Two-photon collisions are the main source of hadron production in the high-energy regime 
of LEP via the process e + e~ — > e + e~7*7* — > e + e~hadrons. Hadrons with high transverse 
momentum are produced by the direct QED process 7*7* — > qq or by QCD processes originating 
from the partonic content of the photon. Next-to-leading order (NLO) QCD calculations are 
available [1,2] for inclusive jet production in quasi-real two-photon interactions. 

The L3 Collaboration published results on inclusive 7r°, Kg [3] and charged hadron [4] pro- 
duction in quasi-real two-photon collisions. The inclusive 7r° and 7r ± differential cross sections, 
measured as a function of transverse momentum, exhibit a clear excess over NLO QCD calcu- 
lations [5] for large transverse momentum. In this Letter, inclusive jet production is studied, in 
similar two-photon interactions, for a centre-of-mass energy of the two interacting photons, W^, 
greater than 5 GeV. The jets are measured in the transverse momentum range 3 < pt < 50 GeV 
and in the pseudo-rapidity interval < 1. The analysis of jet production allows a compar- 
ison of the measurements to NLO QCD predictions, expected to be largely independent of 
fragmentation functions and hadronisation models. 

2 Data and Monte Carlo 

The data used for this analysis were collected by the L3 detector [6] at centre-of-mass energies 
yfs = 189 — 209 GeV, with a luminosity weighted average value of yfs = 198 GeV, and a total 
integrated luminosity of 560 pb -1 . Results on inclusive jet production at LEP for a smaller 
data sample at lower yfs were previously reported [7]. 

The process e + e~ — > e + e~hadrons is modelled with the PYTHIA [8] event generator with 
an event sample two times larger than the data. In this generator, each photon can interact as 
a point-like particle (direct process), as a vector meson (VDM process) or as a resolved photon 
(resolved process), leading to six classes of events. Since both incoming photons are assumed 
to be on the mass shell, PYTHIA is modified to generate the photon flux in the Equivalent 
Photon Approximation [9]. Predictions from the PHOJET [10] Monte Carlo program are 
also compared with the data. The following Monte Carlo generators are used to simulate the 
relevant background processes: KK2f [11] for e + e~— > qq (7); KORALZ [12] for e + e - ^ t + t~(7); 
KORALW [13] for e+e"^ W+W" and DIAG36 [14] for e+e"^ e+e-r+r". Jet hadronisation 
is simulated with the JETSET [8] parton shower algorithm. Events are simulated in the L3 
detector using the GEANT [15] and GHEISHA [16] programs and passed through the same 
reconstruction program as the data. Time dependent detector inefficiencies, as monitored 
during each data taking period, are included in the simulations. 

3 Event selection 

Two-photon interaction events are collected predominantly by the track triggers [17] with a 
low transverse momentum threshold of about 150 MeV. The selection of e + e~ — > e + e~ hadrons 
events [18] consists of: 

• A multiplicity cut. To select hadronic final states, at least six objects must be detected, 
where an object can be a track satisfying minimal quality requirements or a calorimetric 
cluster of energy greater than 100 MeV. 
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• Energy cuts. To suppress background from beam-gas and beam-wall interactions, the 
total energy in the electromagnetic calorimeter is required to be greater than 500 MeV. 
In order to exclude e + e~ annihilation events, the total energy deposited in the calorimeters 
must be less than 0.4 yfs. 

• An anti-tag condition. Events with a cluster in the luminosity monitor, which covers 
the angular region 31 < 9 < 62 mrad, with an electromagnetic shower shape and energy 
greater than 30 GeV are excluded. 

• A mass cut. The mass of all the visible particles of the event, including clusters in the 
luminosity monitor, must be greater than 5 GeV. In this calculation, the pion mass is 
attributed to tracks and electromagnetic clusters are treated as massless. The visible 
mass distribution for data and Monte Carlo is shown, after all cuts, in Figure [T] A wide 
range of masses is accessible. 

About 3 million hadronic events are selected by these criteria. The background level of this 
sample is less than 1% and is mainly due to the e + e~ — > qq(7), e + e~^ t + t~ and e + e~^ 
e + e~T + T~ processes. 

4 Jet definition and composition 

Jets are formed from good quality tracks and electromagnetic clusters. The tracks have a 
transverse momentum greater than 400 MeV, an absolute pseudorapidity less than 1 and a 
distance of closest approach to the primary vertex in the transverse plane less than 4 mm. The 
number of hits must be greater than 80% of the maximum number expected from the track 
angle. For a transverse momentum less than 20 GeV, the momentum and direction of the tracks 
are measured with the central tracker. For the tracks with transverse momentum above 20 GeV, 
the track momenta are replaced with that derived from the energy of their associated cluster 
in the electromagnetic and hadronic calorimeters, assuming the pion mass. Tracks associated 
with muon chamber hits are rejected. An electromagnetic cluster must have an energy greater 
than 100 MeV in at least 2 neighbouring BGO crystals and an absolute pseudorapidity less 
than 3.4. There should be no charged track within an angle of 200 mrad around the cluster 
direction and the associated energy in the hadron calorimeter must be less than 20% of the 
electromagnetic energy. 

Jets are constructed using the k t jet algorithm KTCLUS [19]. This algorithm uses cylindrical 
geometry in which the distance between two objects i, j of transverse momenta pu and ptj is 
defined as dij = min(p^,pf) [(r/i — rjj) 2 + — §j) 2 ]/D 2 where rji and r/j are the pseudorapidities 
of the objects, $j and $j their azimuthal angles with respect to the beam axis and D is a 
parameter of the algorithm which determines the size of the jet. The standard value D = 1 
is used. A distance parameter dk equal to p 2 k is also associated to each object. At the first 
iteration of the algorithm, the objects are the tracks and electromagnetic clusters defined above. 
At each iteration of the algorithm, the d^ and dk are ordered. If the smallest distance is a d^, 
the corresponding objects i and j are replaced by a new object, a "precluster" , formed by 
adding the 4- momenta of the objects i and j. If the smallest distance is a d/~ associated with 
a particle, this is considered as a "beam jet" particle and is removed from the list of objects. 
If the smallest distance is a dk associated with a precluster, this defines a "hard jet" and is 
removed from the list of objects. The procedure is iterated until all objects define beam or 
hard jets. Only hard jets with p t > 3 GeV and < 1 are further considered for this analysis. 
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In Table ^ the data are compared to the Monte Carlo at reconstructed and generated 
levels for: the number of jets, the mean number of jets per event with at least one jet, the mean 
number of particles per jet and outside the jets. For different pt intervals, comparisons are made 
of the mean number of tracks and electromagnetic clusters per jet and of transverse momentum 
of the leading particle divided by that of the jet. The standard deviations of these distributions 
are also quoted. For Monte Carlo at generator level, all particles with mean life time less 
than 3 x 10~ 10 s are allowed to decay and jets are formed from the photons, charged pions, 
charged and neutral kaons, protons and neutrons. Both Monte Carlo programs underestimate 
the number of particles inside and outside the jets. The predicted number of electromagnetic 
clusters is too low for all p t . The amount of energy carried by the most energetic particle of 
the jet is correctly reproduced, except in the highest p t interval. The number of particles per 
jet is shown in Figure El 

Figure El shows the distributions of \rj\ for particles, i.e. clusters and tracks, tracks and jets 
in two intervals of the jet transverse momentum, p t < 20 GeV and p t > 20 GeV. The detector 
acceptance for tracks, calorimetric clusters and jets is well reproduced by Monte Carlo models. 

5 Differential cross section 

The differential cross section for inclusive jet production as a function of p t is measured for 
W 77 > 5 GeV, with a mean value of (W 77 ) — 30 GeV, and a photon virtuality Q 2 < 8 GeV 2 , 
with (Q 2 ) ~ 0.2 GeV 2 . This phase space is defined by Monte Carlo generator-level cuts. 
Results are presented in 9 p t bins between 3 and 50 GeV. 

The p t distribution of the jets is presented in Figure |H The total background is listed 
in Table El Events from the e + e~— > e + e _ r + r~ process dominate the background at low p t 
while hadronic and tau-pair annihilation events dominate it at high p t . To measure the cross 
section, the background is subtracted bin-by-bin. The migration due to the p t resolution is 
corrected by a one-step Bayesian unfolding [20]. The data are corrected for the selection 
efficiency which includes acceptance, and is calculated bin-by-bin as the ratio of the number 
of fully simulated jets selected in PYTHIA over the number of generated jets, as formed by 
the KTCLUS algorithm applied to particles at generator level. The efficiency decreases with p t 
from 61% to 15%. 

The level 1 trigger efficiency is obtained by comparing the number of events accepted by 
the independent track and calorimetric energy triggers [21]. It varies from 97% to 100%. The 
efficiency of higher level triggers is about 98% and is measured using prescaled events. The 
differential cross section and the overall efficiency, which take into account selection and trigger 
efficiencies, are given as a function of p t in Table 121 

Sources of systematic uncertainties on the cross section measurements are the uncertainties 
on the estimation of the selection and trigger efficiencies, the limited Monte Carlo statistics, 
the background subtraction procedure, the selection procedure and the Monte Carlo modelling. 
Their contributions are shown in TableHH The uncertainty due to the selection procedure is eval- 
uated by repeating the analysis with different selection criteria: the multiplicity cut is moved to 
5 and to 7 objects, the requirement on the number of hits of the tracks is moved to 70% of those 
expected, the isolation angle of clusters is moved to 100 mrad, and jets with a particle account- 
ing for more than 90% of the jet transverse momentum are rejected. The sum in quadrature 
of the differences between these and the reference results is assigned as systematic uncertainty 
in Table El Varying other criteria, such as the energy cut, the minimum cluster energy or the 
threshold where the track energy is defined by calorimeters, gives negligible contributions. To 
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evaluate the uncertainty on the Monte Carlo modelling, the selection efficiency is determined 
using only one of the PYTHIA subprocesses: VDM-VDM, direct-direct or resolved-resolved. 
The systematic uncertainty is assigned as the maximum difference between these values and 
the reference Monte Carlo. 

The differential cross sections as a function of \rj\ are uniform within the experimental 
uncertainties for both p t < 20 GeV and p t > 20 GeV, albeit in the latter case these uncertainties 
are large. 

The differential cross section da/dp t is described by a power law function Ap^ B , as expected 
from the onset of hard QCD processes, with B = 3.65 ± 0.07. The result of the fit is shown in 
Figure together with a comparison to Monte Carlo predictions. 

In Figure the data are also compared to analytical NLO QCD predictions [2]. For 
this calculation, the flux of quasi-real photons is obtained using the improved Weizsacker- 
Williams formula [22]. The interacting particles can be point-like photons or partons from the 
7 — > qq process, which evolve into quarks and gluons. The GRV-HO parton density functions 
of Reference 23 are used and all elementary 2^2 and 2^3 processes are considered. The 
parameter A^ 5 ^ is set to 130 MeV. The renormalization and factorisation scales are taken to be 
equal: = M = E t /2[V\. To assign uncertainties, the scale is varied by a factor 1/2 or 2, which 
gives a change in the prediction less than 20%. The results of this calculation agree [2] with 
those described in Reference 24. An additional uncertainty in comparison with NLO QCD, 
which is not considered here, might arise from the modeling of the hadronisation process. In 
a similar study [25] it was evaluated to be below 10% for p t > 10 GeV and decreasing with 
increasing p t . The agreement with the data is poor in the high-pt range, as previously observed 
in the case of inclusive ir° [3] and 7r ± [4] production in similar two-photon reactions. In Figure 
H3 the data are divided in two W 77 ranges, W 7J > 50 GeV and W 17 > 50 GeV and compared 
to the analytical NLO QCD predictions [2]. For W 71 > 50 GeV there is a clear discontinuity 
in the slope near p t = 25 GeV, due to the direct contribution. At high p t , the disagreement 
between data and theoretical calculations is still present. 
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Variable 






Data 




PYTHIA 


PHOJET 












Reconstructed 


Generated 


Reconstructed 


Generated 


Total number of jets 




68792 




107140 


1 88302 


65781 


1 05633 


Number of jets / event 


1.2 


± 0.1 


(0.5) 


1 4 


(0 7) 

V / 


1.3 (0.7) 


1 9 


(0 A) 

V / 


1.1 (0.5) 


N(particles) / jet 
N(particles) outside jets 


6.1 
14.4 


± 0.1 
± 0.1 


(2.5) 
(8.4) 


10.0 


(9 ^ 

(7.0) 


5.3 (2.4) 
13 6 (9 3) 


^ 7 
12.4 


(7 3) 


6.1 (2.4) 
18 4 (8 8) 


N(tracks) / jet 


3<p t <5 GeV 


2.2 


± 0.1 


(1.3) 


2.3 


(1.3) 




2.4 


(1.3) 






5<p t <10 GeV 


2.4 


± 0.1 


(1.3) 


2.6 


(1.3) 




2.8 


(1.4) 






10<p t <25 GeV 


2.5 


± 0.1 


(1.6) 


2.9 


(1.3) 




3.0 


(1.6) 






25<p t <45 GeV 


2.7 


± 0.2 


(1.7) 


3.3 


(1.6) 










N(clusters) / jet 


3<p t <5 GeV 


3.7 


± 0.1 


(2.4) 


2.0 


(2.0) 




3.1 


(2.2) 






5<p t <10 GeV 


3.9 


± 0.1 


(2.6) 


1.8 


(1.9) 




3.3 


(2.4) 






10<p t <25 GeV 


3.9 


± 0.1 


(3.0) 


1.6 


(1.8) 




3.3 


(2.5) 






25<p t <45 GeV 


3.8 


± 0.3 


(3.0) 


1.4 


(1.7) 










p t (leading) / p t 


3<p t <5 GeV 


0.50 ± 0.01 


(0.18) 


0.53 


(0.18) 


0.50 (0.18) 


0.51 


(0.18) 


0.46 (0.17) 




5<p t <10 GeV 


0.54 ± 0.01 


(0.20) 


0.55 


(0.19) 


0.50 (0.20) 


0.52 


(0.19) 


0.43 (0.17) 




10<p t <25 GeV 


0.63 ± 0.01 


(0.23) 


0.60 


(0.20) 


0.48 (0.22) 


0.60 


(0.24) 


0.39 (0.19) 




25<p t <45 GeV 


0.69 ± 0.03 


(0.23) 


0.56 


(0.14) 


0.47 (0.25) 









Table 1: Mean value and standard deviation (in brackets) of multiplicities and p t fractions for the jets in data and Monte Carlo events, 
at generator level as well as after reconstruction. The uncertainties on the mean values are quoted for the data. For Monte Carlo, they 
are always lower than the precision of the last digit. 



Pt 


(Pt) 


Background 


Reconstruction 


Trigger 


da/dp t 






[GeV] 


[GeV] 


[%] 


efficiency [%] 


efficiency [%] 


[pb/GeV] 






3-4 


3.4 


4.6 ± 0.1 


60.8 ± 0.2 


95.8 ± 0.3 


(13 ± 1 ± 1) 


xlO 


i 


4-5 


4.4 


5.6 ± 0.1 


57.2 ± 0.3 


95.9 ± 0.5 


(40 ± 1 ± 3) 






5-7.5 


5.9 


7.8 ± 0.1 


53.2 ± 0.3 


96.2 ± 0.5 


(11 ± 1 ± 1) 






7.5-10 


8.5 


11.1 ± 0.1 


48.9 ± 0.5 


96.6 ± 1.0 


(30 ± 1 ± 2) 


xlO" 


-i 


10-15 


11.9 


14.0 ± 0.2 


44.9 ± 0.6 


96.8 ± 1.4 


(88 ± 3 ± 7) 


xlO" 


-2 


15-20 


17.1 


16.0 ± 0.4 


39.2 ± 0.9 


96.9 ± 2.0 


(30 ± 2 ± 3) 


xlO" 


-2 


20-30 


24.0 


18.6 ± 0.8 


31.6 ± 0.8 


97.3 ± 2.1 


(90 ± 7 ± 8) 


xlO" 


-3 


30-40 


34.1 


18.9 ± 1.5 


20.5 ± 1.3 


97.3 ± 2.5 


(31 ± 5 ± 2) 


xlO" 


-3 


40-50 


44.7 


19.6 ± 1.6 


15.2 ± 1.9 


98.5 ± 2.8 


(11 ± 3 ± 2) 


xlO" 


-3 



Table 2: Background level, reconstruction efficiency, trigger efficiency and differential cross 
section as a function of p t for \rj\ < 1 and W 77 > 5 GeV. The first uncertainty is statistical and 
the second systematic. The average value of p t for each bin, {p t ), is also given. 



Pt 


Trigger 


Monte Carlo 


Background 


Selection 


Monte Carlo 


[GeV] 


efficiency [%] 


statistics [%] 


subtraction [%] 


procedure [%} 


modelling [%] 


3-4 


0.3 


0.3 


< 0.1 


8.4 


0.3 


4-5 


0.5 


0.5 


0.2 


7.0 


1.3 


5-7.5 


0.5 


0.5 


0.3 


6.6 


1.5 


7.5-10 


1.0 


1.0 


0.6 


4.8 


2.4 


10-15 


1.4 


1.3 


0.9 


7.0 


2.6 


15-20 


2.1 


2.4 


1.7 


8.0 


3.3 


20-30 


2.2 


2.6 


2.7 


6.0 


4.8 


30-40 


2.6 


6.4 


5.2 


< 0.1 


6.2 


40-50 


2.8 


12.4 


9.6 


< 0.1 


12.4 



Table 3: Systematic uncertainties on the inclusive jet cross section as a function of p t . 
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Figure 1: Distribution of the visible mass for selected events. The Monte Carlo distributions 
are normalised to the luminosity of the data. Various contributions to the background (back.) 
are shown as cumulative histograms. 
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Figure 2: Distribution of the number of particles per jet for jets with p t > 3 GeV and \r)\ < 
1. The Monte Carlo distributions are normalised to the luminosity of the data. Various 
contributions to the background (back.) are shown as cumulative histograms. 
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Figure 3: Distributions of the pseudo rapidity \rj\ for a) and b) particles and tracks used to 
form jets withp t < 20 GeV andp t > 20 GeV, respectively. "Particles" include both calorimetric 
clusters and tracks, c) and d) distributions of \rj\ for reconstructed jets with p t < 20 GeV and 
p t > 20 GeV, respectively. The Monte Carlo distributions are normalised to the luminosity of 
the data. In a) and b) the higher Monte Carlo lines refer to particles and the lower ones to 
tracks. Various contributions to the background are shown as cumulative histograms in c) and 
d). 
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Figure 4: Distribution of the number of jets with \r]\ < 1 as a function of p t . The Monte 
Carlo distributions are normalised to the luminosity of the data. Various contributions to the 
background (back.) are shown as cumulative histograms. 
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Figure 5: Inclusive jet differential cross section da/dp t a) compared to PYTHIA and PHOJET 
Monte Carlo predictions and the result of a power law fit (solid line); b) compared to NLO 
QCD calculations [2] (solid line). The theoretical scale uncertainty is less than 20%. 
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Figure 6: Inclusive jet differential cross section da/dpt for events with two-photon centre-of- 
mass energy, W 17 , below and above 50 GeV. NLO QCD calculations [2] are superimposed to 
the data. The discontinuity around 25 GeV is due to the direct contribution. 
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